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^ ! Abstract 

We review recent experimental results on D° — D° mixing and CP 
violation charm decays. These studies provide complementary constraints 
on many different extensions of the Standard Model. Observation of CP 
^ ! violation in charm decays at the current level of experimental sensitivity 

would be clear signals of New Physics. 
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1 Introduction 



A process in which a particle changes to its antiparticle (flavor quantum number 
changes by two units) is called particle-antiparticle oscillation or mixing. Particle- 
antiparticle mixing has been observed in all four flavored neutral particle-antiparticle 
systems, i.e., in neutral kaon, both neutral B meson systems and neutral D meson 
system. The particle-antiparticle mixing phenomenon causes an initial (at time t — 0), 
pure D° meson state (below expressions are written for neutral D meson system, 
however they are the same in other three neutral meson systems that exhibit mixing) 
to evolve in time to a linear combination of D° and D° states: 




M 




(1) 



where M and T are 2x2 Hermitian matrices. Diagonal elements of the effective 
Hamiltonian H cff = M — |r describe flavor-conserving transitions D° — > D° and 
D° — >■ D°, while off-diagonal elements describe the flavor- changing transitions D° -h- 
D°. The hermiticity of M and T requires M 12 = M* ± , Ma = M-*, T 12 
Ta = r*j, and the CPT invariance requires M n = M 2 2 = M and r n = r 22 = T 
The eigenstates of the effective Hamiltonian H c g are 



and 



| J D li2 )=p| J D°)±g| J D°), 
while the corresponding eigenvalues are 



(2) 



Ai, 2 = (^m - -rj ± ^ [M 12 - -r 12 j = m 1>2 - -r 1)2 . (3) 

The coefficients p and q are complex coefficients, satisfying \p\ 2 + \q\ 2 = 1, and 
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The real parts of the eigenvalues Ai j2 represent masses, m 12 , and their imaginary 
parts represent the widths r 12 of the two eigenstates |-D 12 ), respectively. The time 
evolution of the eigenstates is given by 



\D 1 , 2 (t))=e 1 , 2 (t)\D 1>2 ), e li2 (t) = e 



■i(mi,2- iri : 2/2)t 



(5) 



The time-dependent decay amplitudes of D° and D° states decaying to final state 
/ are then given by 



(f\n\D°(t)) 
(f\n\W(t)) 



\ \[ei(t) + e 2 (t)\ A f + ^ [ ei (t) - e 2 (t)\ A f j , 
\ [ei(0 - e 2 (t)} A f + [ ei (f) + e 2 (t)} A f j , 



(6) 
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where the instantaneous amplitudes Af and Af are defined as Af = (/|7/|-D°), and 
Af = (/|"H|D°). The time- dependent decay rates are obtained by squaring the above 
decay amplitudes: 

T(D°(t) -> /) = |^/| 2 e- rt f 1 + ' A/|2 cosh(yn) - -Re [A/] sinh^H) 
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— cos(xrt) + Im[Xf] sin(a;rt) 



|2 



2 
9 



1 I I \ 12 

e -rt U /I cosh ( 2/rt ) _ ^ e [ A/ ] S i n h(yrt) 



1 — lAy 



cos(xTt) — Im[Xf] sin(xrt) , (9) 



where A/ = and x and y are dimensionless mixing parameters defined as 

x = l — \ (10) 

Ti — r 2 , . 

v = (ii) 

and T = (ri + T 2 )/2 is the mean decay width. 

In case of a non-zero mass difference (x ^ 0), the mixing is a consequence of a 
pure D° -H- D° transition, and in the case of non-zero lifetime difference (y ^ 0), 
the mixing is a consequence of the shorter-living eigenstate dying out. Since the 
remaining longer-living eigenstate is a linear combination of weak eigenstates D° and 
D°, an initially pure D° sample produces some fraction of D° over time. 

The time evolution of D°, D° — >■ / decays is exponential with the lifetime t d o = 
1/r, modulated by the mixing parameters x and y (see the expressions above). Time- 
dependent measurements of D° and D° decays thus enable us to measure the mixing 
parameters x and y. Since the dependence on x and y is specific for each decay mode, 
different decay modes exhibit different sensitivities to the parameters x and y. 

Out of the four flavored neutral meson systems the neutral D meson system is 
the only one in which down-type quarks are involved in the mixing loop (see Fig. 
[T|). The neutral pion is its own antiparticle and top quark decays before it forms 
a hadron and therefore cannot oscillate. Studies of charm mixing offer therefore a 
unique probe for New Physics (NP) via flavor changing neutral currents in the down- 
type quark sector. In the Standard Model (SM) mixing in neutral D meson system 
can proceed through a double weak boson exchange (short distance contributions) 
represented by box diagrams, or through intermediate states that are accessible to 
both D° and D° (long distance effects), as represented in Fig. HJ Potentially large 
long distance contributions are non-perturbative and therefore difficult to estimate, 
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Figure 1: Short distance (left) and long distance contributions D°-D° mixing in the 
Standard Model. 

so the predictions for the mixing parameters x and y within the SM span several 
orders of magnitude between 10 -8 and 10~ 2 [H |2]. Due to large uncertainties of the 
SM mixing predictions it makes it difficult to identify NP contributions (clear hint 
would be, if x is found to be much larger than y), however, measurements can still 
provide useful and competitive constraints on many NP models, as will be discussed 
later. 

Study of CP violation in decays of charmed hadrons also holds the potential 
for uncovering the NP. In the SM direct CP violation can occur in singly Cabbibo 
suppressed (SCS; c — > dud, c — > sus) decays, but not in Cabbibo favored (CF; 
c —> sud) or doubly Cabbibo suppressed (DCS; c —> dus) decays. This is due to the 
fact that the final state particles in SCS decays contain at least one pair of quark and 
anti-quark of the same flavor, which makes a contribution from penguin-type or box 
amplitudes induced by virtual b-quarks possible in addition to the tree amplitudes. 
However, the contribution of these second order amplitudes are strongly suppressed by 
the small combination of CKM matrix elements V c bV* b . The CP violating asymmetry, 
defined as _ 

a = r(£->/)-r(g->/) 

CP T(D^f) + T(D^f) 1 ; 

is in the SM expected to be at most at the level of 0.1% [3], which is well below 
the current experimental sensitivity. In some NP models the CP asymmetry can be 
significantly enhanced and can be as large as 1% |U El |6j [7] . It is thus widely believed 
that the observation of large CP violation at the order of 1% in charm decays would 
be an unambiguous sign for processes beyond the SM. Direct CP violation occurs 
when the absolute value of the decay amplitude for D to decay to a final state / (A./) 
is different from the one of corresponding CP-conjugated amplitude (Aj). This can 
happen if the decay amplitude can be separated into at least two parts (in case of 
SCS decays the two corresponding SM amplitudes are the tree and the penguin-type) 
associated with different weak and strong phases, Af — \Ai\e tSl e 1 ^ 1 + \A2\e l52 e l< ^ 2 , 
where <pi represents weak phases that switch sign under CP-transformation, and Si 
represent strong phases which are CP-invariant. This ensures that CP-conjugated 
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amplitude, Aj, can differ from Af. Using the definition of Aqp from above one 
obtains Aqp oc sin(4>\ — 4>2)sin{5i — 82). Experimentally, the direct CP- violation is 
probed by measuring the difference between the partial decay widths (T) of D — > f 
and D — > f. In neutral D meson system the time integrated CP asymmetry, Aqp, 
can receive also contributions from CP violation induced by mixing (if \q/p\ 7^ 1) or 
interference between mixing and decay (if arg(qAf/pAf) ^ 0(ir) and / = /). 

2 Recent experimental results on D°-D° mixing 

The most precise constraints on the mixing parameters x and y are obtained using 
the time dependence of D° decays. In time-dependent measurements, the D° decay 
time is calculated as t = m£>o(L ■ Pd°)/\pd°\ 2 > where L is the vector joining the D 0, s 
production and decay vertices, and pp>o and m^o are its momentum and nominal mass. 
Detected tracks of D° decay products are refitted to a common vertex to determine 
the D° decay point, and the production point is taken to be at the primary vertex (in 
pp collisions) or is calculated from the intersection of D° momentum vector with the 
beam spot profile (in e + e~ collisions). Often, the flavor of initially produced neutral 
D mesons needs to be tagged in order to identify D°-D° transitions. The flavor is 
tagged by requiring that neutral D mesons originate from D* + — > D°tt + decays. The 
charge of pion accompanying D° tags the flavor of the neutral D meson at production. 
The energy released in the D* + decay, q = Md*+ — Mpo — m n +, has a narrow peak 
near the threshold for the signal events and thus helps also to reject the background 
candidates. 

Mixing in the D°-D° system has been searched for more than two decades without 
success - until 2007. Three experiments - Belle, BaBar and CDF - so far have found 
evidence for this phenomenon. In the following, the measurements of the mixing 
parameters will be briefly summarized in different D° meson decays. 

2.1 Decays to CP eigenstates 



Belle [8] found first evidence for D°-D° mixing using the ratios of lifetimes extracted 
from a sample of D° mesons produced through the process D* + — > D°tt + , which 
decay to K~tt + , K~K + , or tt~tt + . The time-dependent decay rates of the CF mode, 
K~tt + , and the SCS modes h~h + (h = K or 71") are obtained from time-dependent 
decay rates given in previous section: 



where it has been taken into account that C 1 and |^/^4j| = 1 (|^4j/^4j| <^ 1) 
for D° meson decays to h~h + (K~n + ). The lifetime difference between the CP 



T(D°(t) K-n + ,D°(t) K + n-) oc e 
T{D°{t),D°{t) ^ h + h-) oc e 



*/ T D 

(i+ycp)t/T D0 



(13) 
(14) 
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Figure 2: Proper decay time dependent ratio of D° decays to CP-even eigenstates 
K~K + and 7r~7r + to CP mixed final state K~tt + as measured by Belle [8]. The full 
(red) line shows the dependence of this ratio for measured value of yep, while the 
dashed line shows expectation of no mixing (yep = 0). 



eigenstates h h + and CP-mixed state K n + , yep, is defined as 
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The lifetimes tk-k and r h h are the effective lifetimes extracted from the samples of D° 
mesons decaying to CP mixed final state K~ir + , and CP even final states K~K + 
and 7T~7r + . If \q/p\ = 1 and <fi = 0(ir), the CP symmetry in mixing and interference 
between mixing and decay is conserved, and the parameter y C p corresponds to the 
mixing parameter y. In these time-dependent measurements of neutral D mesons 
decaying to CP eigenstates the indirect CP violation is also probed by comparing 
lifetimes of D° and D° mesons: 



-4, 



t(D°)-t(D°) 
t(D°) + r(D°) 



(16) 



Using sample consisting of around 0.15 (1.2) million reconstructed tagged D° decays 
to h~h + (K~n+) Belle found y CP = (1.13 ± 0.32 ± 0.25)% and A T = (0.01 ± 0.30 ± 
0.15)%. Figure [2] shows proper decay time dependent ratio of D° decays to CP-even 
eigenstates K~K + and 7r~7r + to CP mixed final state K~7t + as measured by Belle 
[8]. Results were later confirmed by BaBar's measurement using D° tagged [9] and 
untagged [TU] samples. 
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2.2 Decays to hadronic wrong sign decays 



Babar found first evidence for D -D mixing in a time dependent study of DCS 
D° — > K + 7i~ decyas [TT]. These decays (also referred to as wrong sign decays) can 
precede both through mixing followed by a CF decay, D° — > D° — > K + tt~ , or directly 
through a DCS decay D° — > K + n~ . To distinguish the two processes, an analysis 
of the decay time distribution is performed. The time-dependent decay rate for the 
two-body wrong sign decays D° — > K + n~ is given by: 



where Rd is a ratio of DCS to CF decays, and the parameters x' and y' are rotated 
mixing parameters, rotated by an unknown strong phase difference between the DCS 
and CF amplitudes, Sk^' x' = x cos 5k* + 2/ sin and y' = ycos5Kn — x sin 5kk- 
The three terms in the time-dependent decay rates of wrong sign decays are due to 
the DCS amplitude, the interference between the DCS and CF amplitudes, and the 
CF amplitude, respectively. BaBar [TT] and CDF [T2] found evidence for oscillations 
in D° — > K + tt~ (samples consisted of around 4 and 13 thousand wrong sign decays, 
respectively) with 3.9 and 3.8 standard deviations, respectively. The most precise 
measurement is from Belle |13], excluding the non-mixing point x' 2 — y' — 0, at 2.1 
standard deviations. Figure |3] shows proper decay time dependent wrong sign ratio 
(R ws = N(D° ->■ K + 7r-)/N(D° ->• K~ir + )) as measured by BaBar p]. Extraction 
of the mixing parameters x and y from the results of these measurements requires 
knowledge of the relative strong phase 5ktt, which can be determined in time inte- 
grated measurements using the correlated D°-D° system produced at ^(3770) [Tlj . 
Both, Belle and BaBar collaborations performed studies of D° and D° samples sepa- 
rately to probe the CP violation in mixing. Results are found to be consistent with 
no CP violation. 

2.3 Three-body decays 

Several intermediate resonances can contribute to a hadronic three-body decay of a 
neutral D meson. For example, D° — > Kg7T + 7r~ decays can proceed via D° — > K*~tt + 
(CF amplitude), D° K° s p° (SCS amplitude and CP eigenstate), D° -> K*+n- 
(DCS amplitude) and many others. In the isobar model, the instantaneous amplitudes 
for D° and D° decays to the three-body final state / are parameterized as a sum of 
Breit-Wigner resonances and a constant non- resonant term (in case of no direct CP 
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Figure 3: The ratio of wrong sign to right sign decays of D° to K + tt~ as measured 
by BaBar [TTJ. The dashed line shows the expected wrong-sign rate as determined 
from the mixing fit. In case of no mixing the ratio is expected to be constant. 

violation, e.g. there is no difference between amplitudes and phases in D° and D° 
decays) : 

A f (s+, s_) = A f (s + , s_) = £ a r e l ^A r (s + , s_) + a NR e^, (19) 

r 

where yfs± is the invariant mass of a pair of final state particles (e.g. Kg7r ), and 
sum runs over all resonances r. The time dependent decay rate for D° decays is thus 
given by (similar expression is obtained for D° decays): 

dT } D ] oc 1^(^,01^^) + |^ a ( a+>a _)| a e-^-»') 

ds + ds-dt 

+2Re[A\(s + , s_)^4.2( s +) s -)} cos (x—\e~r 
+2Im[A 1 (s + , s_)A* 2 (s + , s_)] sin (x-)e~r^ (20) 



T 



(21) 



where Ai^{s + ,sJ) = | LAf(s + ,sJ) ± |^4/(s + ,s_)V Different regions in the s + — s_ 
plane (also called as the Dalitz plot) exhibit different forms of time dependence, as can 
be seen from the above decay rate; therefore, the time-dependent Dalitz plot analysis 
of neutral D meson decays to three-body final state enables us to measure the x and 
y parameters simultaneously. In case the analysis is performed separately for D° and 
D° samples CP violation can be probed by measuring q/p (amplitude and phase) 
directly. This method of measuring mixing parameter x and y was pioneered by Cleo 
in D° — >■ Ks7t + tt~ decays [15], and applied also by Belle in D° — > K^-k + -k~ decays 
[T6] and BaBar in D° -> K^tt + tt~ and D° ->■ K° S K + K~ decays |17]. Measurements 
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Figure 4: Central value (point) and C.L. contours (including statistical, systematic 
and amplitude model uncertainties) in the x — y plane for C.L. = 68.3%, 95.4%, 99.7% 
as obtained by BaBar in [17]. The no-mixing point is shown as a plus sign (+) and 
is disfavored at 1.9 standard deviations. 

are consistent with each other and together provide the most accurate determina- 
tions of x and y (average of all three measurements assuming no CP violation, as 
obtained in Ref. [T8]): x = (0.42 ±0.21)% and y = (0.46 ± 0.19)%. Belle's (BaBar's) 
measurement disfavor the no-mixing hypothesis with a C.L. equivalent to 2.2 (1.9) 
standard deviations. Belle also performed search for CP violation in mixing in their 
measurement and found \q/p\ = (0.86±0.30±0.09), arg(q/p) = (-0.24±0.30±0.09), 
to be consistent with 1 and 0, respectively. Figure H] shows the confidence level (C.L.) 
contours in the x — y plane as obtained by BaBar collaboration. 

Large fractions of D° — > K$K + K~ decays proceed via D° — > Kg(f) (CP-odd) and 
D° — > Kga (980) (CP-even) decays. Belle took advantage of this fact and performed 
Dalitz plot integrated measurement of yep mixing parameter using untagged sample 
of D° — y K$K + K~ decays (19]. They measure the effective lifetimes of D° mesons, 
Ti on, OFF? in two different regions of K + K~ invariant mass (at peak (ON) and </> 
sidebands (OFF)), which are given by t n,off = (1 + (1 — 2foN,OFF)ycp) T D°, where 
/on, off is CP-even fraction in ON or OFF region. The obtained value of yep is 
consistent with y obtained by Babar's Dalitz plot analysis (T7J using these decays 
only. The sensitivity to y of both approaches is similar. 

BaBar found evidence of D°-D° mixing also with time-dependent analysis of wrong 
sign D° K + 7T~7r decays [20]. Analysis is similar to analysis of wrong sign 
K + 7i~ decays, however since the strong phase difference, Sktttt , varies across the 
available three-body phase space the Dalitz plot analysis needed to be performed. 
The time dependent decay rate depends on the DCS amplitude Aj and CF amplitude 
Aj and is given by: 



dT(s + , s ,t) 
ds + ds dt 



{\Aj(s + , s ) | 2 + \ Aj(s + , s )Aj(s + , s )\(ycosSj - xcos6j)t/r(22) 
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Figure 5: (a) Proper time distribution for right sign events with the fit result super- 
imposed. The distribution of background events is shown by the shaded histogram, 
(b) Proper time distribution for wrong sign events, (c, d) sq and s + projections with 
superimposed fit results (line). 

+\^ 7 (s + ,s )\ 23 ^(t/r)^e^, 

where 5j(s + ,s ) = arg\Aj(s + , so)Aj(s + , s )\. Figure [5] shows projections of proper 
decay time and Dalitz plot variables s and s + with superimposed fit results. The 
no-mixing point is excluded with a significance of 3.2 standard deviations. 

2.4 World average and constraints on New Physics models 

All existing measurements of D°-D° mixing (for full list of existing measurements 
check [18]) are combined by the Heavy Flavor Averaging Group by performing a 
global x 2 fit. The average central values are found to be x = (0.631q;2o)% an d 
y = (0.75±0.12)% [18]. Figure |6] shows the resulting 1 — 5 standard deviation contours 
in x — y and \q/p\ — arg(q/p) planes. No D°-D° mixing hypothesis is excluded at 
10.2 standard deviations, while there is no hint for CP violation at current level of 
sensitivity. 

Golowich et al. [21] studied implications of existing D°-D° measurements on 
many NP models. They find in many scenarios strong constraints that surpass those 
from other search techniques and provide an important test of flavor changing neutral 
currents in the up-quark sector. One simple extension to the SM that they studied 
is the addition of a fourth family of fermions. The obtained constraint on the CKM 
mixing parameters V c yV* b , (b' is the down-quark of the fourth generation) is an order 
of magnitude stronger than those obtained from unitarity considerations of the CKM 
matrix. 
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Figure 6: Two dimensional contours for mixing parameters x and y (left) and CP 
violating parameters \q/p\ and arg(q/p) (right) as obtained by HFAG when averaging 
all relevant existing measurements. 

3 Recent experimental results on direct CP viola- 
tion in charm 

In this section we will focus only on time-integrated measurements searching for direct 
CP violation in decays of charmed hadrons since previous section already covered re- 
sults from time-dependent studies. Searches of CP violation were performed in past 
15 years in over 30 decay modes of D°, D + , and D$ mesons by Belle, BaBar, Cleo, 
CDF, FOCUS, E796 and E687 experiments (the full list of all direct CP violation 
measurements is available in Ref. p3]). No evidence for CP violation were found so 
far, however the measurements have only started reaching interesting level of sensi- 
tivity below 1% in some decay modes. In order to increase the sensitivity to or below 
0.1% level not only larger samples but also very good control over the systematic 
uncertainties will be needed. These uncertainties are dominated by the uncertainties 
in asymmetries in the detection and reconstruction of particles of opposite charge. 
In addition, forward-backward production asymmetries, Afb, resulting from Z° and 
virtual photon interference and higher order loops in the production of c-anti-c quark 
pairs in e + e~ collisions, results in asymmetries in the distribution of D decay prod- 
ucts in regions of varying efficiency in the detector. Estimation of these factors used 
to rely upon Monte Carlo simulated studies, with questionable assumptions about 
charge dependent interaction effects, resulting in systematic uncertainties in Acp's 
in the 1 — 5% range. In the past years, new insights in using real data rather than 
simulations have led to reduction of these uncertainties to the 0.2 — 0.5% range. These 
corrections and systematic uncertainties are decay mode dependent, however it is im- 
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portant to note that are determined by the statistics of the samples used and will 
thus decrease with increasing sample sizes. 

3.1 Direct CP violation in neutral D meson decays to CP 
eigenstates 

Time-integrated CP asymmetry, Acp, defined in Eq. [12] receives in case of neutral 
D meson decays contributions of direct and indirect (mixing induced) CP violation. 
The latter is independent of the decay mode (universal) and its time-dependence is 
for decays to CP eigenstate given by [22J: 

A cpif) ~ \ncp (y {\p/q\ ~ cos<p + x (\p/q\ + \q/p\) sin<p) -. (23) 

Integrating over time yields 

A CP = 4p + / A CP (t)D(t)dt « af P + ^a*£. (24) 

The factor — in front of the second term, indirect CP violation term, is 1, if the 
reconstruction efficiency does not depend on proper decay time, D{t) = 1 (e.g. as 
it is the case at Belle and Babar). At CDF and LHCb the displaced-track trigger 
requirements reject candidates with short decay times which results in — > 1. 

Recently, CDF presented measurements of Aqp in neutral D meson decays to CP 
eigenstates K~K + and 7r~7r + [23J. The D° mesons are required to originate from 
D* + D°7T + decays. The raw reconstructed asymmetry, defined as A veco (hh) = 
(N(D° ->■ hh) - N(D° -> hh))/(N(D° ->■ hh) + N(D° ->■ hh)), is a sum possible 
CP asymmetry in D° decays and detector induced slow pion reconstruction efficiency 
asymmetry: 

D* + -> D°7if -> [h + h-]7ij => A veco {hh)* = A CP {hh) + A £ (tx s ). (25) 

In order to correct for the latter asymmetry and to extract the Acp(hh) asymmetry 
CDF reconstructed also tagged and untagged samples of D° — > K~tt + decays. The 
raw reconstructed asymmetries in these two samples are given by: 

D* + -)• D°tt+ -)• [K-tt + ]tt+ A reco {K7T)* = Acp(Ktt) + A e {ir s ) + A £ (Kix){2Q) 
D°^K-n + => A reco {Kn) = A C p{Ktt) + A £ {Ktt).{27) 

(28) 

The asymmetry Acp(hh) is then obtained by the following combination of raw re- 
constructed asymmetries: 

A CP {hh) = A ICCO {hh)* - A vcco {Ktc)* + A vcco (Ktx), (29) 
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Figure 7: Comparison of Acp(7nr) (left) and Acp(KK) (right) in the parameter space 
{ci 1 CPi a 'cp) measured by CDF, Belle and Babar collaborations. 

assuming that the production asymmetry to be (which is the case since at Tevatron 
the initial state, pp, is charge symmetric), the efficiency to reconstruct D* can be 
factorized into tt s and D° reconstruction efficiencies and that kinematic equations are 
equal across all four samples. The CDF collaboration found, 

A CP (nny = (+0.22 ±0.24 ±0.11)%, (30) 
A C p(KK)* = (-0.24 ±0.22 ±0.10)%, (31) 

which are consistent with CP conservation. These are the most precise measure- 
ments up to date of time-integrated asymmetries of D° decays to K~K + and tt~tt + 
final states. Figure [7] shows the comparison of this measurement with the ones per- 
formed by Belle and Babar collaborations [2U [25] in the parameter space (a™p,a^p) 
(according to Eq. |24"|) . 

Recently the LHCb collaboration presented a measurement of the difference of 
the time-integrated CP asymmetries in D° — > K~K + and D° — > 7r~7r + decays us- 
ing data collected in 2010 [26] in which production asymmetry and detector related 
asymmetries cancel: 

A Teco (KK)* = A CP (KKy + A £ (n s ) + A P (D*), (32) 
A rcco (7T7r)* = A CP (<K7cy + M<K s ) + Ap(D*), (33) 
AA CP = A CP (KK)* - A C p(tttt)* = A rcco (KK)* — A rcco (7r7r)*. (34) 

They measured AA C p = (—0.28 ± 0.70 ± 0.25)% to be consistent with 0. Using also 
tagged and untagged samples of D° — > K~tt + decays the production asymmetry of 
D° mesons, Ap(D°), can be expressed and measured if world averages for Acp(hh) 
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Figure 8: Measured production asymmetry of D° mesons in bins of rj (left) and pt 
(right). 

are taken as external inputs. Figure El shows the measured D* production asymmetry 
as a function of transverse momentum, px, and pseudo-rapidity, rj. At current level 
of statistics the production asymmetry does not significantly depend on pt and rj, 
but when averaged out over LHCb's acceptance region is inconsistent with 0. 

3.2 Direct CP violation in decays of charged D mesons 

Recently, Belle collaboration presented results of the measurement of Aqp asymmetry 
difference between SCS D + — > 07r + and CF — > (f)n + decays in the region of 
\M(K + K~) - rn^l < 16 MeV/c 2 [27J. In case of charged D {s) decays to 0tt+ the 
reconstructed asymmetry can be divided into four contributions 

A reco = A CP + A FB (cos6*) + MK + K~) + A £ ( Pnj cosO*), (35) 

where A FB (cos9*) is the forward-backward production asymmetry of the mesons, 
A £ (K + K~) is the asymmetry in the reconstruction of K + K~ pair, and A e (p n , cos9 v ) 
is the asymmetry in the pion detection. Naively, one would expect the A e (K + K~) 
to be zero, since the K + K~ pair is charge symmetric. This would be the case, if 
the K + and K~ momentum spectra would be identical. As can be seen in Fig. 
the momentum spectra are different which leads to non-zero A £ (K + K~). This 
asymmetry can be very precisely determined using A £ (K) (measured on data with 
— > (j)7i + and D° — > K~tt + decays) and normalized kaon phase space distributions 
P(x) (x = (p, cos6))\ 

A £ (K + K-) = J(P K+ (x)-P K ~(x))A £ (K)dx, (36) 
AA £ (K + K-) = A°(K + K~) -A°°(K + K~) = (+0.111 ±0.025)%. (37) 

The difference in reconstructed asymmetries is according to Eq. [35] then given by 
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Figure 9: Background subtracted momentum distributions of same and opposite 
charge relative to that of the D (a) or D s (b) decays. 

AA TCCO = A CP {D + + AAfb{cos9*) + AA £ (K + K~), (38) 

where the intrinsic Acp in CF Z)+ — > <j)7r + is assumed to be negligible and provided 
the measurement is performed in bins of the (cos9*,p n , cosO^) phase space. To obtain 
Acp{D + — > 07r + ) and AAfb{cos9*) the weighted average of Ay4 reco in each bin of 
cos9* is calculated and then A C p(D + — > (f)n + ) and AA FB (cos6*) are extracted by 
adding or subtracting bins at ±cos8*. The CP asymmetry in D^<p7T + decays is 
found to be (0.51 ± 0.28 ± 0.05)%, and is consistent with no direct CP violation in 
these decays. This is the most precises measurement of Aqp in these decays up to 
date. In addition no significant difference between forward-backward asymmetries in 
the production of the D + and mesons is found as can be seen in Fig. [TOj 

4 Conclusions 

Studies of charm mixing and CP violation offer unique opportunities to search for 
processes beyond the SM. Large samples of neutral D meson collected by Belle, BaBar 
and CDF experiments enabled us to find evidence for D meson mixing after 31 years 
of its discovery. Next round of experiments (LHCb at CERN, Bellell at KEK in 
Japan, and SuperB in Italy) will collect even orders of magnitude larger samples 
which might reveal also CP violation in charm. The future was never more bright 
for charm physics. 
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